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Abstract 

We consider a tiling model of the two-dimensional square-lattice, where each site is tiled with 
one of the sixteen Wang tiles. The ground states of this model are all quasi-periodic. The systems 
undergoes a disorder to quasi-periodicity phase transition at finite temperature. Introducing a 
proper order-parameter, we study the system at criticality, and extract the critical exponents 
characterizing the transition. The exponents obtained are consistent with hyper-scaling. 
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It has been known for a while that the two- or three-dimensional space may be tiled by 
ordered but aperiodic tilings, in addition to periodic lattice structures This aperiodic 
order is realized in nature by certain alloys, called quasicrystals {2!, 3|, which are believed to 
exhibit in thermal equilibrium aperiodic crystalline order [J, - long-range positional order 
lacking any periodicity. Aperiodic tiling models are extensively used to analyze quasycrystals 
at zero or finite temperature. However, understanding of the transition region between the 
disordered (fluid) phase and the quai crystal phase is incomplete yet. 

Here we consider a model of interacting tiles on a square, two-dimensional, lattice. The 
model has been previously studied in jg, 7|. Each site of the square lattice is tiled with a tile, 
characterized by four labels attached to its edges. The labels take one of possible six labels 
(or colors). The interaction is with nearest neighbor tiles, and the bond energy is zero if the 
labels of both neighboring edges match, or one otherwise. It was found by Ammannjs] that if 
one limits the allowed tiles to a group of 16 tiles (out of all the possible 6*^ tiles), then all zero 
energy states of the model (also known as perfect tilings) are non-periodic. The 16 Amman 
tiles are presented in figure [H The perfect tilings are then the ground states of this system. 
Their non-periodicity can be shown using a mapping of the six tile labels into one of two 
symbols: 5* for labels {1,2} and L for labels {3,4,5,6}. Since all Amman tiles have labels 
belonging to the same set {S or L) on both horizontal (vertical) edges, the mapping thus 
classifies the 16 tiles into 4 types according to their set along each axis: {1} S-S, {2, 3, 4, 5} 
S-L, {6, 7, 8, 9} L-S and {10, 11, 12, 13, 14, 15, 16} L-L. Identifying the symbols S and L with 
the short (S) and long (L) Fibonacci tiles, it follows from the properties of the Amman tiles 
that any perfect tiling is mapped into a two-dimensional square Fibonacci tiling {q], thus 
aperiodic. 

n 

The finite temperature behavior of this tiling model was studied numerically in |6|- The 
model has multiple ground states (uncountable infinite number for the infinite plane), all 
are aperiodic, and thus its dynamics upon fast cooling was suggested to be a model for 
glassiness. The lattices studied in jslwere in the range 8 < < 32 (A^ being the lattice 
linear size) with free boundary conditions. Numerical results supported the existence of 
a phase transition, measured by a growing peak in the specific heat, and the transition 
was concluded to be of second order. Recently Q], it was shown that the phase transition 
observed in [6| is a disorder (fiuid) to quasicrystal transition. Phase transition analysis in 
[3] followed an analytical approach supported by numerical simulations. The transition was 
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studied using the overlap of a configuration with a ground state 7. The fraction of tiles in 
a configuration c matching a ground state 7 is denoted 0(c, 7). The overlap, normalized by 
its averaged over all ground states ipic) = J (f>{c, 'j)dX{'j), is then thermally averaged to yield 

^^^^^ = zJ — W) ' 

where the integration is over configurations c (ensemble averaging), H{c) is the energy of the 
configuration, and Z is the partition function Z = J exp(— (c)). The high temperature 
limit of Q 13 is, by normalization, 1. 

To account for the infinite number of ground states, Koch and Radin defined the quantity 
Qrk 

qRKW) = -J Qpil) HQp{^)]dX{^), (2) 

and then qRK^P « (3c) = 0. Analytical calculation showed that qrk vanishes identically 
for sufficiently high finite temperatures. Numerical simulations were then used to show qjiK 



does not vanish 
persented in 



1 for low temperatures, thus proving the existence of a transition. Simulations 
7| employed fixed boundary conditions corresponding to a specific ground 
state a. Then, the overlap of the configuration with the chosen ground state a was used 
to approximate qnK at low temperatures, where it is expected to contribute dominantly. 
Looking at system sizes 32 < < 256, they concluded that the transition is of third or 
higher order. In addition, it was suggested that the transition has no renormalization fixed 
point. 

Square Fibonacci quasicrystals are known to exhibit distinct delta-function peaks in the 
fourier transform, in a similar fashion to the peaks observed in crystalline solid j^. We 
therefore suggest to use the weights of these distinct peaks to define an order parameter for 
the disorder-to-quasicrystal transition, one that is simpler and easier to access numerically 
then qRK- Many frequencies show a peak, all are irrational and related to the 'golden mean' 
r = (1 + \/5)/2. Frequency chosen for this study is ko = ((r — l)/r, (r — 1)/t), but results 
are similar for all related frequencies. The (complex) order parameter is then the amplitude 
of the peak at /cq 

q = I e^'o %,irfr. (3) 

Tile 1 is chosen for simplicity, as it is the only tile which is the only one of its type, i.e. it 
is the only tile with both edges being S edges. Similar results are obtained using any other 
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FIG. 1: (Color Online) The sixteen Amman tiles. The four rows correspond to the four types of 
tiles (see text). 

type e.g tiles 2-5, 6-9 or 10-16. 

In the following we present numerical simulations for 10 < < 400 and free boundary 
conditions. Simulations starts from a ground state and then thermalized at the desired 
temperature. Thermalization during cooling was found to be significantly less efficient. As 
expected, the order-parameter auto-correlation times measured near criticality are extremely 
long, and thus we took each run to last 10^ Monte-Carlo Steps. To account for multiple 
ground states and to assure good cover of phase space we repeated the simulations starting 
from different ground states. Finite size scaling analysis was used to calculate the critical 
temperature and critical exponents of the phase transition. 

In order to determine the critical temperature, we analyzed the Binder cumulant 



Q = l- 



(3|gP) 



where (■ ■ ■) denotes ensemble average. Crossings of Q for system sizes 10 < < 100, 
presented in figure [21 show no significant finite size effects within our accuracy. The relatively 
large error estimates due to the long relaxation processes allow for a moderate accuracy in 
fixing the critical temperature, which is estimated to be = 0.4216(5) (here and below 
the number in paranthesis is the uncertainty in the last digit). The critical temperature 
bund is in agreement with previous estimates: Tc = 0.42(1) Q| and /3c = 1/Tc ~ 2.4 
3]- We then studied larger lattices at two temperatures near criticality, Ti = 0.4211 and 
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FIG. 2: (Color Online) Binder cumulant of the tiling model. The collapse of crossings for various 
lattice sizes signatures the location of the critical point. 

T2 = 0.4222. Finite size scaling of thermodynamic quantities at Tc provides estimates for 
the critical exponents. In the following we verify that the critical exponents obtained for Ti 
and T2 are similar, in order to ensure that our inaccuracy in does not take us out of the 
critical regime for the lattice sizes studied. The order parameter |g| is expected to scale as 
{\q\){L) ~ L^^/^. Figure [3] presents numerical results from measurements at Ti and T2 and 
fits to the power-law form. Based on the two fits we estimate fi/v = 0.30(3). Similarly, the 
susceptibility defined by 



X 



N\{W)-{\q\?) 
T 



(4) 



was also measured at Ti and T2 and fitted to the scaling form x ~ -^^^'^ (figure HI). Esti- 
mation of 7/z/ from the two fits leads to 7/z/ = 1.40(5). Note that these two independent 
measurements satisfy the hyperscaling relation 2/3/z/ + 7/z/ = d. 

In order to estimate the critical exponent we studied the derivative of the binder cumu- 
lant {Q) with respect to the inverse temperature /3. Based on finite-size-scaling arguments 

dQ 



Q' 



d(3 



T=T, 



is expected to diverge like L^^'^. It is easy to see that this derivative is obtained from the 
energy and order-parameter moments in the following form 
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FIG. 3: (Color Online) Order parameter measurements performed at Ti and T2. Power law behavior 
is fitted, L=10 is excluded from the analysis. Fitted exponents are 0.28 and 0.33. 




FIG. 4: (Color Online) Susceptibility x fitted to a power law behavior at Ti and T2. Fitted 
exponents are 1.44 and 1.36. 

Measurements presented in figure [5] indeed show a power law behavior of Q'{N), and the 
exponent estimated is 1/z/ = 1.05(15). This value for u is consistent with measurements of 
the specific heat 



(5) 

i.e. a = 2 — du = 0, or u = 1. The 



which are best-fitted by a logarithmic growth (figure [6 
results are consistent with the analysis performed in (6|, leading to v = 1.6(5). 

In conclusion, we show that the use of the Bragg peak allows for an analysis of the disorder 
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FIG. 5: (Color Online) Measurements of Q' for Ti = 0.4211. Fitted exponent is 0.99. 
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FIG. 6: (Color Online) Specific heat measurements at T\ and T^- Data is well fitted by a logarithmic 
function. For comparison, best power law fits are presented in the inset, yielding exponents 0.23 
and 0.25. 

to quasicrystal transition in a two-dimensional lattice model based on Amman tiles. The 
transition occurs at a finite temperature and is of a second order. Critical exponents were 
measured and shown to satisfy hyperscaling relations. This model is therefore suitable for 
study of the critical emergence of quasi-periodic order. 
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